Pionic susceptibility for charged pions in asymmetric nuclei 
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At low energies the particle-hole {NN~^) part of the pionic susceptibility in isospin-symmetric 
nuclear matter is known to behave very differently from the susceptibility in finite nuclei due to the 
presence of an energy gap in the NN~^ excitation spectrum. In this note we show that for charged 
pions in N ^ Z nuclei the changes due to the gap are very similar to those in the symmetric case, 
except at very low momenta, where a qualitatively different behavior is found. 
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1 In several applications [1, 2] the polarization function 
of a nuclear medium is needed for very low energy trans- 
fer w ~ 0. As an example, we will consider here the pionic 
susceptibility xi^jk) [3], but our discussion applies like- 
wise to other isovector channels. For finite nuclei, x(^i k) 
is often calculated in a local-density approximation using 
results from the homogenous medium. In this case the 
particle-hole (NN^^) contribution, which is dominant at 
low energies, is described by the Lindhard function. How- 
ever, when applied to finite nuclei, this approximation is 
justified only if the energy uj is much larger than the 
lowest-lying iVA^~^ excitation energy. The modifications 
due to this energy gap in the excitation spectrum were 
discussed in Ref. [1] for the isospin-symmetric case and 
are summarized below. In this note we will focus on the 
isospin-asymmetric case. 

For charged pions in isospin-asymmetric matter the 
original Lindhard function [4, 5] cannot be used. Instead, 
it is convenient to split the Lindhard function into direct 
and crossed contributions. We define 
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where a and b indicate neutron (n) or proton (p), p^p"^ 
and p^p^ denote the corresponding Fermi momenta, and 
m is the nucleon mass. For completeness, the explicit 
expression for (t>ab-^ is given in the appendix. Using this 
definition, e.g., the non- interacting (Fermi-gas) suscepti- 
bility for a 77+ is given by 
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X^^l (lo, k) = 2^ [</.p„-l {u, k) + 4>np-^ {-UJ, k)] , (2) 

where / denotes the ttNN coupling constant. The sus- 
ceptibility for a 7r~ can be obtained by replacing co <-> 
—UJ. In realistic cases one has to include corrections 
from the short-range NN correlations (Ericson-Ericson 
Lorentz-Lorenz effect [3]) which lead to an RPA resum- 

mation of x^r^ ■ 

In symmetric matter the indices a and b can be omit- 
ted. In this case it is obvious from Eq. (1) that, for uj 0, 



the function (f> fulfills (j){uj 7^ 0, fc 0) = as a conse- 
quence of Pauli blocking, whereas for lu = one finds 
(j){LJ — 0, k ^ 0) = — mpi?/(27r^) . The reason for this 
behavior is that in the latter case both numerator and 
denominator of Eq. (1) vanish. However, as pointed out 
in Ref. [1], this results in a bad approximation for the 
susceptibility of a finite nucleus. Due to the finite exci- 
tation energy of the lowest-lying NN^^ excitation, the 
denominator cannot vanish for uj ^ 0. As a consequence 
the NN~^ part of the susceptibility of a finite nucleus 
must go to zero as fc — > 0. Following Ref. [1], this ef- 
fect is easily included by adding an energy gap A to the 
NN~^ excitation energy (fc^ -I- 2k ■ p)/(2m). Then the 
pionic susceptibility can be written as 
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A, fc) -f (^(-cj - A, fc)] . (3) 



The drastic effect of this modification can be seen from 
the dashed lines in Fig. 1 which are in agreement with 
Fig. 2 of Ref. [1] . The short-dashed line corresponds to 
the usual Lindhard function, Eq. (2), whereas the long- 
dashed line represents the Lindhard function with a gap, 
Eq. (3). 

We now return to the asymmetric case. To be spe- 
cific, we will discuss the case of neutron-rich matter. 
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In this case 4>pn-i{uj,k) is non-zero 



for all values of uj and k, since Pauli blocking can never 
be complete. This is in contrast to (/)„p-i (w, fc), which 
is identically zero due to Pauli blocking for all mo- 
It should also be noted that 



menta k < p^p"^ 
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4>pn-i [uJ — 0, fc) and thus also x^^l {'^ = 0, k) are complex 
for k < 2p'p\ 

The previous discussion is valid for infinite matter 
which is adequate e.g., for the interior of a neutron star. 
In view of the strong modifications in the symmetric case 
due to the energy gap in finite nuclei, it is interesting to 
include similar corrections also in the asymmetric case. 
We will consider two modifications: 

1. A neutron star is stabilized against /3 decay by the 



Fermi energy of the electrons, e 
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FIG. 1: NN ^ part of the pionic susceptibility at nuclear 

matter density (q,i + Qp ~ 0.17 fm~^) for fixed = as func- 
tion of k. Dashed lines; symmetric matter, corresponding to 
Pf^ = Pf' ~ 1.36 fm~^. Short dashes: usual Lindhard func- 
tion, Eq. (2). Long dashes: modified Lindhard function with 
gap in the NN~^ excitation spectrum, Eq. (3). Dash-dotted 
and solid lines: asymmetric matter with Qn = l-^Qp, cor- 
responding to p^' = 1.45 fm^^ and p^' = 1.26 fm^^. Dash- 
dotted line: usual Liudliard function, generalized to the asym- 
metric case, Eq. (3). Solid line: with gap in the excitation 
spectrum and proton energies shifted by Vc, adjusted such 

that e(?> = e<r^ 



In fact, over a wide range of momenta it seems to be a 
very accurate approximation to neglect the asymmetry 
completely. However, at momenta below » — p'p 
the susceptibility stays almost constant in the asymmet- 
ric case (solid line), while it goes to zero for A; ^ in the 
symmetric case (long-dashed line). As mentioned above, 
this qualitative difference stems from the fact that in the 
asymmetric case Pauli blocking can never be complete. 
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APPENDIX: EXPLICIT EXPRESSIONS 

In the explicit expression for the function (f)g^),-i{ui,k) 
defined in Eq. (1) four different cases must be distin- 
guished: 

(a) k > p^p^ +P^p (no Pauli blocking): 

</'afe-(c^, k) = 4^1, {oj, k, -pf,pP) , (A.l) 



The Coulomb potential is zero as a consequence of 
charge neutrality. This situation is very different in 
finite nuclei, where the whole Fermi sea of the pro- 
tons is shifted upwards with respect to the Fermi 
sea of the neutrons by the Coulomb potential Vc, 
such that for /3-stable nuclei the Fermi energies of 
protons and neutrons are almost equal (see, e.g.. 
Fig. 1(a) of Ref. [6]). Therefore we shift all proton 

energies by Vc ~ {p^p^ ^ — Pp^ ^)/(2m) . 

2. As in the symmetric case, we include the energy 
gap A in order to account for the discrete NN~^ 
excitation spectrum at low energies. It should be 
noted, however, that even if there is a gap in the 
pp~^ or nn~^ excitation spectrum, there is no gap 
in the pn~^ spectrum of /3-unstable nuclei. 

Including both effects, the pionic susceptibility for a 7r+ 
reads 
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+ <P„p-i{-u-A + Vc,k)]. (4) 

The result of these modifications concerns the imaginary 
part of XTr+i^jk), which now vanishes for ui = 0. The 
real part changes in a way very similar to the symmet- 
ric case, as can be seen from Fig. 1, where the dashed- 
dotted line corresponds to the susceptibility of infinite 
asymmetric matter, Eq. (2), while the solid line contains 
the energy gap A and the Coulomb potential Vc, Eq. (4). 



(b) \pP-pP\<k<pP+pP: 

+ f^U^,k,p^^'> ~k,p^^), (A.2) 

{c) k<\pP -Poland pP <pf: 

(/)„6-i (w, k) =/^"li iuJ,k, -pP, -pP - k) 

+ filU^,k,-pP-k,pP-k) 

+ fi;U^,k,pP ~k,pP), (A.3) 

(d) k < \pp^ — pp^ I and pp^ > pp^ (complete Pauli 
blocking): 

0„6-i(a;,fc)=O. (A.4) 

In these equations fjj^-i and denote the integrals 

over p|| and p± in the regions of py where the integration 
over p± is unrestricted or restricted by Pauli blocking, 
respectively: 

A^) ( 1 . m 

fa,-d^,k,p„p,) = ^^ 

X [2k{pi -P2){k'^ - k{pi +P2) - 2mw) 

+ {Ak^pP ' - (fc' - 2muj^f) Inx] , (A.5) 
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+ ipP^ -p'-p^^ + 2rmj)lnx\ , 

with 

+ 2kpi — 2mui — is 
k'^ + 2kp2 — 2mco — is 



If one is only interested in the real part of Eq. (2), the 
Pauli-blocking effects in the two terms cancel, and it is 
(A. 6) possible to use the simple formula (A.l) for all cases, see 
Eqs. (10.6) to (10.8) of Ref. [7]. However, for Eqs. (3) 
and (4), the distinction of the four cases (a) to (d) is 
necessary. 

(A.7) 
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